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Abstract 

Isaac  Newton  proposed  hanging  a  bucket  of  water  by  a  cord  in  the 
Principia.  If  the  cord  is  twisted  and  the  bucket  is  then  released,  it 
begins  to  spin  and  the  surface  of  the  water  acquires  a  paraboloidal 
shape.  In  this  paper,  the  parabolic  profile  as  a  function  of  the  angle  of 
rotation  is  derived,  as  well  as  the  period  of  the  torsional  oscillations  as 
a  function  of  the  initial  parameters  of  the  system. 

1.  Introduction 

If  a  vessel,  hung  by  a  long  cord,  is  so  often  turned  about  that  the  cord  is 
strongly  twisted,  then  filled  with  water,  and  held  at  rest  together  with  the 
water...  [then]  while  the  cord  is  untwisting  itself...  the  vessel,  by 
gradually  communicating  its  motion  to  the  water,  will  make  it  begin 
sensibly  to  revolve,  and  recede  by  little  and  little  from  the  middle,  and 
ascend  to  the  sides  of  the  vessel,  forming  itself  into  a  concave  figure  (as  I 
have  experienced)  and  the  swifter  the  motion  becomes,  the  higher  will  the 
water  rise  ...  This  ascent  of  the  water  shows  its  endeavor  to  recede  from 
the  axis  of  its  motion;  and  the  true  and  absolute  circular  motion  of  the 
water  ...  becomes  known,  and  may  be  measured  by  this  endeavor.  [1] 

Newton’s  bucket  is  well  known  to  philosophers  of  science,  who 
have  pondered  the  metaphysics  of  why  the  liquid  in  a  rotating  container 
adopts  a  curved  surface.  Ernst  Mach,  for  example,  postulated  that  the 
parabolic  shape  must  be  due  to  the  existence  of  matter  in  the  universe, 
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specifically  of  the  distant  “fixed  stars”  relative  to  which  the  bucket  rotates 
[2],  However,  the  simple  physical  questions  are:  How  does  the  shape  of 
the  water  surface  vary  as  the  cord  supporting  the  bucket  unwinds?  What  is 
the  period  of  the  torsional  oscillations? 

2.  Surface  Profile  From  a  Force  Analysis 

Suppose  water  partly  fills  a  right  cylindrical  bucket  (of  radius  R ) 
that  is  rotating  about  its  vertical  z  axis  of  symmetry  at  angular  speed  (O. 
Make  the  key  assumption  (to  be  discussed  later)  that  the  liquid 
instantaneously  follows  the  motion  of  the  bucket.  (The  speed  of  the  bucket 
is  restricted  to  be  less  than  some  value  comax  that  prevents  both  the 
spinning  water  from  spilling  over  the  top  edge,  and  the  bottom  of  the 
bucket  from  being  exposed  at  the  axis  of  rotation.)  Denote  the  cylindrical 
coordinates  of  any  point  on  the  surface  of  the  water  as  (r,  (/),  z )  where  the 
origin  lies  on  the  axis  at  the  bottom  of  the  bucket.  The  axes  are  fixed  in  the 
laboratory  frame  and  do  not  rotate  with  the  bucket.  The  angular  speed  of 
the  water  and  bucket  is  co  =  d(j)  /  dt .  Two  forces  act  on  a  bit  of  water  at  the 
surface.  One  is  gravity  vertically  downward.  The  other  can  be  alternatively 
described  as  being  due  to  the  pressure  from  the  surrounding  water  [3],  as  a 
buoyant  force  [4],  or  simply  as  a  normal  force  [5],  The  resultant  of  the  two 
forces  is  a  radially  inward  centripetal  force  [6]  in  the  inertial  laboratory 
frame,  or  equivalently  a  radially  outward  centrifugal  force  [7]  in  the 
noninertial  frame  of  the  bucket.  By  considering  the  vertical  and  horizontal 
components  of  Newton’s  second  law  [8],  one  finds  that  the  water  surface 
adopts  the  paraboloidal  shape 

2  2 

CO  r  ... 

z  =  zo  +  ~ —  (!) 

2g 

where  g  =  9.80  m/s  is  Earth’s  gravitational  field  strength  and  z0  is  the 
height  of  the  water  at  the  center  of  the  spinning  bucket.  (Another  way  to 
derive  this  result  is  to  note  that  the  surface  of  the  water  must  be  an 
equipotential  relative  to  the  sum  of  the  gravitational  and  centrifugal 
potential  energies  [9].)  The  value  of  z0  can  be  related  to  the  total  mass  m 
of  water  in  the  bucket, 

R 

m  =  2np  J  rzdr  (2) 

o 
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where  p  =  1000  kg/m  is  the  density  of  water.  Substituting  Eq.  (1)  into  (2) 
and  performing  the  integral  gives 


m  -  npR" 


z0+- 


co2R2^ 


4  g 


(3) 


Denote  the  height  of  the  water  in  the  bucket  when  it  is  stationary  by  h. 
Then  putting  (0=  0  in  Eq.  (3)  implies 


m  =  rcpR  h, 


(4) 


so  that  Eq.  (3)  can  be  rearranged  as 


zo 


co2R2 
4  g 


(5) 


Substituting  this  result  into  Eq.  (1)  leads  to  the  normalized  expression 


*=1- 


co 


co¬ 


ir 


2  \ 


R- 


(6) 


where  comax  =  2R~l  ^[gh  .  Note  that  z  =  0  at  r  =  0  when  co  =  comm ,  in 
agreement  with  the  parenthetical  discussion  of  rumax  above  Eq.  (1).  Also 
note  that  zth  =  2  at  r  —  R  when  co-  fttniax ,  which  implies  that  the  bucket 
must  initially  be  fdled  no  more  than  halfway  with  water,  to  prevent  liquid 
from  spilling  out  at  the  maximum  angular  speed.  Equation  (6)  is  plotted  in 
Figure  1  for  three  different  values  of  co!  (Omm.  For  any  angular  speed, 
z  =  h  when  r  /  R  =  2  .  As  the  bucket  spins  faster,  the  water  level  drops 

in  the  center  and  rises  up  near  the  walls,  as  Newton  noted. 
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Figure  1.  Profiles  of  the  water  surface  for  three  different  angular  speeds. 


3.  Angular  Speed  From  an  Energy  Analysis 

The  moment  of  inertia  of  the  water  in  the  spinning  bucket  is 


Jr2 2 nprzdr  ■ 


1  +  - 


co 


3co" 


(V) 


using  Eqs.  (4)  and  (6),  where  the  moment  of  inertia  of  the  water  when  the 
bucket  is  at  rest  is  70  =  mR  /  2  .  The  moment  of  inertia  increases  as  water 
is  flung  farther  away  from  the  axis  of  rotation  with  increasing  angular 
speed,  up  to  a  maximum  value  of  4/0  /  3 . 

Just  as  the  elastic  potential  energy  of  a  spring  with  a  particle 
attached  to  its  end  is  he  /  2  where  x  is  the  translational  displacement  of 
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the  particle  and  k  is  the  spring  constant,  so  the  torsional  potential  energy  of 
the  rope  with  the  bucket  attached  to  its  end  is 

UT=\c(f)2  (8) 


where  (f)  is  the  angular  displacement  of  the  bucket  and  c  is  the  torsional 
constant  of  the  rope.  The  gravitational  potential  energy  of  the  water 
relative  to  the  bottom  of  the  bucket  is 


R  z 


R  ( 
2 


UG  =  J  J  2  nprdzdrgz  =  xpgj  rh 


o  o 


co 2  larr2  ^ 

1 - ^  + 


v 


,  yL  2  n2 

^rnax  ^rnax  ^ 


dr  (9) 


using  Eq.  (6).  With  the  help  of  Eq.  (4),  the  integral  simplifies  to 


U, 


mgh 


G 


1  +  - 


03 


.4  h 


3  co 


(10) 


max  j 


which  reduces  to  the  expected  result  if  co-  0 .  Note  the  seemingly 
paradoxical  fact  that  even  though  the  water  is  cylindrically  symmetric  and 
UG  is  therefore  independent  of  (j)  as  measured  in  the  rotating  frame  of  the 
bucket,  UG  is  a  function  of  co  which  in  turn  depends  on  the  angle  (j)  as 
measured  in  the  inertial  frame  of  the  laboratory.  The  resolution  of  this 
paradox  is  that  the  angular  acceleration  is  presumed  to  be  small  enough 
that  z  can  be  taken  to  be  independent  of  (f)  over  any  2n  range,  and  yet  the 
height  of  the  water  at  a  given  radius  varies  over  the  course  of  many 
revolutions  of  the  bucket.  As  noted  in  Ref.  [7],  water  in  a  9-cm-diameter 
Lucite  cylinder  spinning  at  a  constant  rate  of  300  rpm  takes  about  1 
minute  to  attain  its  equilibrium  paraboloidal  shape,  indicating  that  the 
coupling  between  z  and  (f>  is  weak  but  nonzero. 

Assume  that  the  mass  of  the  bucket  is  negligible  compared  to  that 
of  the  water.  Then  the  total  potential  energy  U  of  the  system  is  the  sum  of 
UT  and  UG.  Suppose  that  the  rope  is  twisted  through  an  initial  angle  </>0  and 
the  bucket  is  released  from  rest,  so  that  the  initial  kinetic  energy  is  Kx  =  0 
and  the  initial  potential  energy  is 

Ui=\c$  +\mgh.  (11) 


When  the  rope  has  untwisted  to  some  angle  (f)  so  that  the  bucket  is  rotating 
at  angular  speed  co,  the  kinetic  energy  is 
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vu/mav 


(12) 


according  to  Eq.  (7),  and  the  potential  energy  is 


1 


1 


Uf=-c(p  +-mgh  + 


mgh 

6co* 


co  . 


(13) 


Conservation  of  energy  now  implies  that 

c($-<l>2)  =  Il>a?  +  Io‘^\+mgha>4. 


(14) 


Noting  that  IQcqaax=2mgh,  we  can  rearrange  Eq.  (14)  into  the 
normalized  form 


4  2 

CD  .  CO 
- (-2 — — 


cfi b 


co„ 


CO, 


mgh 


(15) 


v  '  u  / 

Solving  this  biquadratic  equation  gives 
,2 


or 


or. 


l  +  £ 


-1 


(16) 


where 


P=  C^° 


mgh 


(17) 


The  dimensionless  constant  (5  is  the  ratio  of  the  initial  torsional  potential 
energy  ccp{]  /  2  to  the  initial  gravitational  potential  energy  mgA  /  2 .  The 
square  root  of  Eq.  (16)  is  plotted  in  Figure  2  for  three  different  values  of  / 3 . 
Note  that  the  maximum  value  of  (5  is  3  if  co  is  not  to  exceed  &>max  when  the 
cord  has  fully  unwound  at  0  =  0.  Furthermore,  even  at  the  midpoint  of  the 
bucket’s  oscillations  when  0  =  0,  increasing  /3  from  1  to  3  increases  co  by 
only  55%. 
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Figure  2.  Angular  speed  of  the  water  as  a  function  of  the  fractional  unwinding  of  the 
rope  for  three  different  initial  numbers  of  twists  of  the  bucket. 


4.  Period  of  Torsional  Oscillations 


Take  the  square  root  of  Eq.  (16),  substitute  co  =  d(j)l  dt,  and 
separate  variables.  Then  integrate  over  a  quarter  period  77 4  as  the  bucket 
passes  through  its  equilibrium  position  and  the  cord  fully  winds  back  up, 
to  get 


1 


i +p 


fa  J 


-1/2 


77  4 


dtp  £ymax  J  dt . 


(18) 


Multiply  both  sides  of  this  equation  by  p  .  Then  make  the  change  of 
variable  in  the  left-hand  integral  to  0  where  sin  0  =  (f>  /  (fo .  Perform  the 
right-hand  integral  and  substitute  Eq.  (17)  into  it  to  eliminate  [3.  Using 
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&>max  =  ^2 mgh  / 10  ,  Eq.  (18)  gives  the  period  of  the  bucket’s  torsional 
oscillation  as 


dO. 


(19) 


The  square  root  in  the  denominator  inside  the  square  brackets  is 
approximately  1  +  0. 5/? cos-  6  in  the  limit  as  [5  — >  0.  Denoting  the  period 
as  T0  in  this  small-angular-amplitude  limit,  one  immediately  obtains 


(20) 


as  expected,  since  /  — >  /0  for  small  co  according  to  Eq.  (7).  The  integral  in 
Eq.  (19)  can  be  numerically  evaluated  for  nonzero  /3,  but  it  is  found  to 
only  increase  slowly  with  (5.  Even  at  the  maximum  value  of  J3  =  3 ,  the 
period  of  oscillation  is  merely  12%  larger  than  T0. 

In  any  case,  Eq.  (16)  gives  an  exact  solution  in  phase  space, 
whereby  quantities  are  expressed  in  terms  of  the  twist  angle  0  rather  than 
in  terms  of  the  elapsed  time  t.  For  example,  substituting  Eq.  (16)  into  (6) 
gives  the  height  of  the  water  at  any  point  in  the  bucket  as  a  function  of  the 
angle  that  the  cord  has  unwound.  In  particular,  at  the  walls  of  the  bucket 
where  r  =  R ,  let  Z  denote  the  height  of  the  water.  Then  the  fractional  rise 
in  the  height  of  the  water  at  the  walls  above  the  stationary  level  is 


Z-h 

h 


l  +  P 


1- 


<t> 


2  A 


-1 


(21) 


which  is  equal  to  the  normalized  square  of  the  angular  speed  of  the  bucket, 
according  to  Eq.  (16).  As  already  mentioned,  Z  =  2h  when  J3  =  3  and 
0  =  0. 


5.  Closing  Remarks 

Why  is  the  period  12%  longer  for  large-angle  oscillations  of  this 
torsional  pendulum  than  it  is  for  small  amplitudes?  The  reason  is  not  the 
same  as  for  a  simple  pendulum.  For  a  simple  pendulum,  the  period 
increases  because  the  approximation  sin  0  ~  0  breaks  down  at  large 
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angles.  Instead,  the  reason  here  is  the  increase  in  the  moment  of  inertia  of 
the  water,  in  accordance  with  Eq.  (7).  In  particular,  if  we  froze  the  water, 
then  the  period  would  be  independent  of  amplitude,  just  as  it  is  for  a  mass 
on  a  Hookean  spring. 

Finally,  let’s  return  to  the  key  assumption  underlying  the  analysis. 
The  viscosity  of  the  water  must  be  high  at  the  walls  and  bottom  of  the 
bucket  if  the  fluid  is  to  instantaneously  adjust  to  the  motion  of  the  solid 
container.  At  the  same  time,  the  viscosity  needs  to  be  low  within  the  bulk 
of  the  fluid  to  prevent  differences  in  angular  speed  between  one  region  of 
the  water  and  another.  Fortunately,  simulations  for  the  spin  of  an 
incompressible  fluid  in  a  cylindrical  container  suggest  that  there  are  viscid 
boundary  layers  in  the  water  near  the  solid  surfaces  of  the  cylinder, 
accompanied  by  an  inner  inviscid  core  [10].  The  situation  is  similar  to 
laminar  flow  over  an  airplane  wing,  with  drag  motion  close  to  the  wing 
and  potential  flow  far  away  from  it. 
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